A hierarchical Bayesian model for inference of copy number variants and their association to gene expression

A number of statistical models have been successfully developed for the analysis of high-throughput data from a single source, but few methods are available for integrating data from different sources. Here we focus on integrating gene expression levels with comparative genomic hybridization (CGH) array measurements collected on the same subjects. We specify a measurement error model that relates the gene expression levels to latent copy number states which, in turn, are related to the observed surrogate CGH measurements via a hidden Markov model. We employ selection priors that exploit the dependencies across adjacent copy number states and investigate MCMC stochastic search techniques for posterior inference. Our approach results in a unified modeling framework for simultaneously inferring copy number variants (CNV) and identifying their significant associations with mRNA transcripts abundance. We show performance on simulated data and illustrate an application to data from a genomic study on human cancer cell lines.Comment: Published in at http://dx.doi.org/10.1214/13-AOAS705 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

The Sodium Channel B4-Subunits are Dysregulated in Temporal Lobe Epilepsy Drug-Resistant Patients

Temporal lobe epilepsy (TLE) is the most common type of partial epilepsy referred for surgery due to antiepileptic drug (AED) resistance. A common molecular target for many of these drugs is the voltage-gated sodium channel (VGSC). The VGSC consists of four domains of pore-forming α-subunits and two auxiliary β-subunits, several of which have been well studied in epileptic conditions. However, despite the β4-subunits’ role having been reported in some neurological conditions, there is little research investigating its potential significance in epilepsy. Therefore, the purpose of this work was to assess the role of SCN4β in epilepsy by using a combination of molecular and bioinformatics approaches. We first demonstrated that there was a reduction in the relative expression of SCN4B in the drug-resistant TLE patients compared to non-epileptic control specimens, both at the mRNA and protein levels. By analyzing a co-expression network in the neighborhood of SCN4B we then discovered a linkage between the expression of this gene and K+ channels activated by Ca2+, or K+ two-pore domain channels. Our approach also inferred several potential effector functions linked to variation in the expression of SCN4B. These observations support the hypothesis that SCN4B is a key factor in AED-resistant TLE, which could help direct both the drug selection of TLE treatments and the development of future AED

Nano-Electrochemical Characterization of a 3D Bioprinted Cervical Tumor Model

Current cancer research is limited by the availability of reliable in vivo and in vitro models that are able to reproduce the fundamental hallmarks of cancer. Animal experimentation is of paramount importance in the progress of research, but it is becoming more evident that it has several limitations due to the numerous differences between animal tissues and real, in vivo human tissues. 3D bioprinting techniques have become an attractive tool for many basic and applied research fields. Concerning cancer, this technology has enabled the development of three-dimensional in vitro tumor models that recreate the characteristics of real tissues and look extremely promising for studying cancer cell biology. As 3D bioprinting is a relatively recently developed technique, there is still a lack of characterization of the chemical cellular microenvironment of 3D bioprinted constructs. In this work, we fabricated a cervical tumor model obtained by 3D bioprinting of HeLa cells in an alginate-based matrix. Characterization of the spheroid population obtained as a function of culturing time was performed by phase-contrast and confocal fluorescence microscopies. Scanning electrochemical microscopy and platinum nanoelectrodes were employed to characterize oxygen concentrations - a fundamental characteristic of the cellular microenvironment - with a high spatial resolution within the 3D bioprinted cervical tumor model; we also demonstrated that the diffusion of a molecular model of drugs in the 3D bioprinted construct, in which the spheroids were embedded, could be measured quantitatively over time using scanning electrochemical microscop

A Systems Biology Approach Reveals a Calcium-Dependent Mechanism for Basal Toxicity in Daphnia magna.

This document is the Accepted Manuscript version of a Published Work that appeared in final form in Environmental Science & Technology, copyright © American Chemical Society after peer review and technical editing by the publisher.The expanding diversity and ever increasing amounts of man-made chemicals discharged to the environment pose largely unknown hazards to ecosystem and human health. The concept of adverse outcome pathways (AOPs) emerged as a comprehensive framework for risk assessment. However, the limited mechanistic information available for most chemicals and a lack of biological pathway annotation in many species represent significant challenges to effective implementation of this approach. Here, a systems level, multistep modeling strategy demonstrates how to integrate information on chemical structure with mechanistic insight from genomic studies, and phenotypic effects to define a putative adverse outcome pathway. Results indicated that transcriptional changes indicative of intracellular calcium mobilization were significantly overrepresented in Daphnia magna (DM) exposed to sublethal doses of presumed narcotic chemicals with log Kow ≥ 1.8. Treatment of DM with a calcium ATPase pump inhibitor substantially recapitulated the common transcriptional changes. We hypothesize that calcium mobilization is a potential key molecular initiating event in DM basal (narcosis) toxicity. Heart beat rate analysis and metabolome analysis indicated sublethal effects consistent with perturbations of calcium preceding overt acute toxicity. Together, the results indicate that altered calcium homeostasis may be a key early event in basal toxicity or narcosis induced by lipophilic compounds

A computational framework for gene regulatory network inference that combines multiple methods and datasets

<p>Abstract</p> <p>Background</p> <p>Reverse engineering in systems biology entails inference of gene regulatory networks from observational data. This data typically include gene expression measurements of wild type and mutant cells in response to a given stimulus. It has been shown that when more than one type of experiment is used in the network inference process the accuracy is higher. Therefore the development of generally applicable and effective methodologies that embed multiple sources of information in a single computational framework is a worthwhile objective.</p> <p>Results</p> <p>This paper presents a new method for network inference, which uses multi-objective optimisation (MOO) to integrate multiple inference methods and experiments. We illustrate the potential of the methodology by combining ODE and correlation-based network inference procedures as well as time course and gene inactivation experiments. Here we show that our methodology is effective for a wide spectrum of data sets and method integration strategies.</p> <p>Conclusions</p> <p>The approach we present in this paper is flexible and can be used in any scenario that benefits from integration of multiple sources of information and modelling procedures in the inference process. Moreover, the application of this method to two case studies representative of bacteria and vertebrate systems has shown potential in identifying key regulators of important biological processes.</p

Gene signatures in wound tissue as evidenced by molecular profiling in the chick embryo model

<p>Abstract</p> <p>Background</p> <p>Modern functional genomic approaches may help to better understand the molecular events involved in tissue morphogenesis and to identify molecular signatures and pathways. We have recently applied transcriptomic profiling to evidence molecular signatures in the development of the normal chicken chorioallantoic membrane (CAM) and in tumor engrafted on the CAM. We have now extended our studies by performing a transcriptome analysis in the "wound model" of the chicken CAM, which is another relevant model of tissue morphogenesis.</p> <p>Results</p> <p>To induce granulation tissue (GT) formation, we performed wounding of the chicken CAM and compared gene expression to normal CAM at the same stage of development. Matched control samples from the same individual were used. We observed a total of 282 genes up-regulated and 44 genes down-regulated assuming a false-discovery rate at 5% and a fold change > 2. Furthermore, bioinformatics analysis lead to the identification of several categories that are associated to organismal injury, tissue morphology, cellular movement, inflammatory disease, development and immune system. Endothelial cell data filtering leads to the identification of several new genes with an endothelial cell signature.</p> <p>Conclusions</p> <p>The chick chorioallantoic wound model allows the identification of gene signatures and pathways involved in GT formation and neoangiogenesis. This may constitute a fertile ground for further studies.</p

Systems biology reveals how altered TGF beta signalling with age reduces protection against pro-inflammatory stimuli

<div><p>Osteoarthritis (OA) is a degenerative condition caused by dysregulation of multiple molecular signalling pathways. Such dysregulation results in damage to cartilage, a smooth and protective tissue that enables low friction articulation of synovial joints. Matrix metalloproteinases (MMPs), especially MMP-13, are key enzymes in the cleavage of type II collagen which is a vital component for cartilage integrity. Transforming growth factor beta (TGFβ) can protect against pro-inflammatory cytokine-mediated MMP expression. With age there is a change in the ratio of two TGFβ type I receptors (Alk1/Alk5), a shift that results in TGFβ losing its protective role in cartilage homeostasis. Instead, TGFβ promotes cartilage degradation which correlates with the spontaneous development of OA in murine models. However, the mechanism by which TGFβ protects against pro-inflammatory responses and how this changes with age has not been extensively studied. As TGFβ signalling is complex, we used systems biology to combine experimental and computational outputs to examine how the system changes with age. Experiments showed that the repressive effect of TGFβ on chondrocytes treated with a pro-inflammatory stimulus required Alk5. Computational modelling revealed two independent mechanisms were needed to explain the crosstalk between TGFβ and pro-inflammatory signalling pathways. A novel meta-analysis of microarray data from OA patient tissue was used to create a Cytoscape network representative of human OA and revealed the importance of inflammation. Combining the modelled genes with the microarray network provided a global overview into the crosstalk between the different signalling pathways involved in OA development. Our results provide further insights into the mechanisms that cause TGFβ signalling to change from a protective to a detrimental pathway in cartilage with ageing. Moreover, such a systems biology approach may enable restoration of the protective role of TGFβ as a potential therapy to prevent age-related loss of cartilage and the development of OA.</p></div

A systems biology approach sheds new light on Escherichia coli acid resistance

In order to develop an infection, diarrhogenic Escherichia coli has to pass through the stomach, where the pH can be as low as 1. Mechanisms that enable E. coli to survive in low pH are thus potentially relevant for pathogenicity. Four acid response systems involved in reducing the concentration of intracellular protons have been identified so far. However, it is still unclear to what extent the regulation of other important cellular functions may be required for survival in acid conditions. Here, we have combined molecular and phenotypic analysis of wild-type and mutant strains with computational network inference to identify molecular pathways underlying E. coli response to mild and strong acid conditions. The interpretative model we have developed led to the hypothesis that a complex transcriptional programme, dependent on the two-component system regulator OmpR and involving a switch between aerobic and anaerobic metabolism, may be key for survival. Experimental validation has shown that the OmpR is responsible for controlling a sizeable component of the transcriptional programme to acid exposure. Moreover, we found that a ΔompR strain was unable to mount any transcriptional response to acid exposure and had one of the strongest acid sensitive phenotype observed

MMP13 and TIMP1 are functional markers for two different potential modes of action by mesenchymal stem/stromal cells when treating osteoarthritis

Mesenchymal stem cells (MSCs) have been investigated as a potential injectable therapy for the treatment of knee osteoarthritis, with some evidence of success in preliminary human trials. However, optimization and scale‐up of this therapeutic approach depends on the identification of functional markers that are linked to their mechanism of action. One possible mechanism is through their chondrogenic differentiation and direct role in neo‐cartilage synthesis. Alternatively, they could remain undifferentiated and act through the release of trophic factors that stimulate endogenous repair processes within the joint. Here, we show that extensive in vitro aging of bone marrow‐derived human MSCs leads to loss of chondrogenesis but no reduction in trophic repair, thereby separating out the two modes of action. By integrating transcriptomic and proteomic data using Ingenuity Pathway Analysis, we found that reduced chondrogenesis with passage is linked to downregulation of the FOXM1 signaling pathway while maintenance of trophic repair is linked to CXCL12. In an attempt at developing functional markers of MSC potency, we identified loss of mRNA expression for MMP13 as correlating with loss of chondrogenic potential of MSCs and continued secretion of high levels of TIMP1 protein as correlating with the maintenance of trophic repair capacity. Since an allogeneic injectable osteoar therapy would require extensive cell expansion in vitro, we conclude that early passage MMP13+, TIMP1‐secretinghigh MSCs should be used for autologous OA therapies designed to act through engraftment and chondrogenesis, while later passage MMP13−, TIMP1‐secretinghigh MSCs could be exploited for allogeneic OA therapies designed to act through trophic repair